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Abstract 

We consider an equation modeling the evolution of a viscous liquid thin 
film wetting a horizontal solid substrate destabilized by an electric field 
normal to the substrate. The effects of the electric field are modeled by 
a lower order non-local term. We introduce the good functional analysis 
framework to study this equation on a bounded domain and prove the 
existence of weak solutions defined globally in time for general initial data 
(with finite energy). 
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1 Introduction 

In this paper, we construct solutions for a thin film type equation with a destabi- 
lizing singular integral term. This term models the effects of an electric field (see 
[19]). From the analytical point of view, this paper belongs to the large body 
of literature devoted to the thin film equation with destabilizing terms such as 
long-wave unstable thin film problems [71 [51 [50] or the Kuramato-Sivashinsky 
equation in combustion and solidification [121 111] . 

More precisely, we are considering the following equation, which is intro- 
duced by Tscluiko and Papageorgiou in jTOj (sec also pTBl): 

Ut + {u^{cu^^-au-XI{u)),^)^^Q xe[0,L], t>0 (1) 

(in |19| , ([T]) is supplemented with periodic boundary conditions) . This equation 
models the evolution of a liquid thin film (of height u) wetting a horizontal solid 
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substrate which is subject to a gravity field and an electric field normal to the 
substrate. The term XI{u) models the effects of the electric field on the thin 
film. The operator I{u) is a nonlocal elliptic operator of order 1 which will be 
defined precisely later on (for now, we can think of it as being the half-Laplace 
operator: /(it) = — (— A)-^/^u). When A > 0, it has a destabilizing effect (it has 
the "wrong" sign). The term au accounts for the effects of gravity, and it is 
also destabilizing when a < ("hanging film"). In [19], it is proved that despite 
these destabilizing terms, positive smooth solutions of ([1]) do not blow up and 
remain bounded in for all time. 

As mentioned above, there are many papers devoted to the study of thin film 
equations with destabilizing terms. In its simplest form, the thin film equation 
reads: 

Uf + {.f{u)Uxxx)x = 0, (2) 

The existence of non-negative weak solutions for ([2]) was first established by 
F. Bernis and A. Friedman [3] for j(u) = m", n > 1. Further results (existence 
for n > and further regularity results) were later obtained, by similar technics, 
in particular by E. Beretta, M. Bertsch and R. Dal Passo [2] and A. Bcrtozzi 
and M. Pugh [5J|5]. Results in higher dimension were obtained in particular by 
Griin in [Ml [131 ^M. ■ The thin film equation with lower order destabilizing terms 
has also received a lot of interest. In particular in [7l|8], the following equation 

+ (/(")"2;i^2; - 5('^)'^x)i^ = (3) 

is considered. Such a destabilizing term (which, unlike that of ([T|), is a local 
term of order 2) models, for instance, the effects of gravity for a hanging thin 
film, or van der Waals type interactions with the solid substrate. In [7|, the 
nonlinearities /(u) = m" and g(u) = u™ are considered and it is proved (among 
other things) that there is no blow-up for m < n -I- 2. In [8|, for /(u) = u and 
g(u) = w"', it is proved that there is blow-up for to > 3 and initial data in iJ^(R) 
with negative "energy" . The reader is referred to [8] for a precise statement. 

In our equation ([TJ, the nonlinearities in front of the stabilizing and desta- 
bilizing terms are the same (/(m) ~ 9iu) = u^), but the destabilizing term is 
elliptic of order 3 and is nonlocal in space. It is known (see [19]) that positive 
smooth solutions of ([ij do not blow up. The goal of the present paper is to 
prove the existence of global in time weak solutions for ([1]) . 

Note that besides the existence of solutions, many important properties of 
the thin film equation ([2]) have been investigated (finite speed of propagation of 
the support, waiting time phenomenon, existence of source- type solutions etc.). 
The key tools in many of these studies are various delicate integral inequalities 
(in particular the so-called a-entropy inequalities and local entropy and energy 
inequalities. See [H [31 [HI [7] ) • It is not clear that similar functional inequalities 
holds for ([1}. One reason is that the algebra involving the operator I{u) is 
considerably more difficult than that of the Laplace operator. Another reason, 
is the obvious difficulty in deriving local estimates (due to the nonlocal nature 
of the operator I{u)). For that reason, we only address the existence issue in 
this paper. 
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As in [3 [S] , the main difficulty in proving the existence of solutions for ([T]) 
comes from the fact that the energy (see ([8]) below) can take negative values. 
In order to obtain a priori estimate, one thus has to use the conservation of 
mass which, for non-negative solutions, gives a global in time bound for the 
solution (see Lemma [T]) . 

With such an estimate in hand, the existence of global in time solutions 
should follow from the construction of approximated solutions satisfying the 
right functional inequalities. Typically, one needs to regularize the mobility 
coefficient u^. One way to proceed is to replace the coefiicient with +£ so 
that the equation becomes strictly parabolic. However, for such a regularized 
equation one cannot show the existence of non-negative solutions (the maximum 
principle does not hold for fourth order parabolic equations) and Lemma [T] is 
of no use. An alternative regularization is to replace the mobility coefficient 
with a function fe{u) which satisfies in particular fe{u) ~ je. For such (more 
degenerate) mobility coefficient, solutions are expected to be strictly positive 
and therefore smooth. This second regularization procedure was first suggested 
by Bernis and Friedman [3] and is used in particular by Bertozzi and Pugh [51[7]. 
However, the local in time existence for such a degenerate equation is not clear 
to us, since the corresponding proofs in [J] rely on Schauder estimates which are 
not classical (and perhaps tedious) with our non-local singular term l(u). For 
this reason, we choose the first regularization approach; this requires us to pay 
attention to the lack of positivity of the approximated solutions. In particular, 
the norm is not controlled, and the norm will be controlled by combining 
the energy inequality with the entropy inequality. The idea of combining the 
energy together with the entropy in order to get a Lyapunov functional appeared 
previously in |17| where a thin film equation with a nonlinear drift term is 
thouroughly studied in all dimensions. 

The main contribution of this paper is thus to introduce the precise func- 
tional analysis framework to be used to treat the term l(u) and to provide a 
method for constructing weak solutions satisfying the proper a priori estimates. 

Rather than working in the periodic setting, we will consider equation ([1]) 
on a bounded domain with Neumann boundary conditions (these Neumann 
conditions can be interpreted as the usual contact angle conditions and seem 
physically more relevant - the periodic framework could be treated as well with 
minor modifications). Further details about the derivation of ([1]) will be given 
in Section [5] Since the gravity term is of lower order than the electric field term, 
it is of limited interest in the mathematical theory developed in this paper. We 
will thus take 

a = and c = A = 1. 
We thus consider the following problem: 

wt + (/(u)(uxx - /(w))x)^ = fori-eO, t>0 

= 0, /(w)(w,,-/(u)), =0 for a; ear!, t>0 (4) 
u(a;,0) = uq(x) for a; S 57. 
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The domain is a bounded interval in M; in the sequel, we will always take 
= (0, 1). The mobility coefficient /(u) is a function / : [0, +oo) (0, +oo) 
satisfying 

f{u) - u" as li ^ (5) 

for some n > 1. The operator / is a non-local elliptic operator of order 1 
which will be defined precisely in Section [3] as the square root of the Laplace 
operator with Neumann boundary conditions (we have to be very careful with 
the definition of / in a bounded domain) . 



A priori estimates. As for the thin film equation ([2]), we prove the existence 
of solutions for ^ using a regularization/stability argument. The main tools 
are integral inequalities which provide the necessary compactness. Besides the 
conservation of mass, we will see that the solution m of Q satisfies two important 
integral inequalities: We define the energy £{u) and the entropy e{u) by 

= ^ / {ul{t) +u{t)Iu{t))dx and e{u){t) = ^ / G{u{t))dx 

where G is a non- negative convex function such that fG" = 1. Classical solu- 
tions of ^ then satisfy: 

E(u){t)+ [ f fiu)[{u,,-I{u))x]^dxds<£{uo), (6) 
Jo Jn 

e{u){t) + J J {uxx)'^ dx ds + J J Uxl{u)x dx ds < e{uo). (7) 

Similar inequalities hold for the thin film equation ([2])- However, we see here 
the destabilizing effect of the nonlocal term I{u): First, we note that as in [71 H] 
the energy £{u) can be written as the difference of two non- negative quantities: 

£{u){t)^\\u{t)\\l,^^^^-\\uit)\\l,^^^, (8) 

and may thus take negative values. Similarly, the entropy dissipation can be 
written as 

{uxxfdxds+l I Uxl{u)xdxds = \\u{t)\\\^. - \\u{t)\\'^i , 

so the entropy may not be decreasing. 

Nevertheless, it is reasonable to expect ^ to have solutions that exist for 
all times. Indeed, as shown in [H], the conservation of mass, the inequality ^ 
and the following functional inequality (see Lemma [T} 

\\u\\l,^^^<a£{u)+l3\\u\\l,^^^, \/ueH\n), 

implies that non-negative solution remains bounded in L°°{0,T; H^{^1)) for all 
time T. 
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Furthermore, the interpolation inequahty 

\\u{t)\\ls < C\\uit)\\fj,\\uit)\\^, 



2" '"H- 2 
yields 

If* , . , . 1 



<n){t) + - \\u{r)r^,dr<e{ua) + - \\u{r)r^, dr, 
and so the entropy remains bounded for all time as well. 

Main results. We now state the two main results proved in this paper. They 
should be compared with Theorems 3.1 and 4.2 in [3J pp.l85&194]. The first 
one deals with non-negative initial data whose entropies are finite. 
We recall that G is a non-negative convex function such that 

= -r^ for all u > 0. 

Theorem 1. Let n > 1 and uq E H^(n) be such that uq > and 

G{uo) dx < oo. (9) 

o 

For all T > there exists a function u(t, x) > with 

u e C{0,T; L^n)) n L°°{0,T; H\n)), e L^[Q,T-Hl{n)) 
such that, for all <j) £ 2?([0, T) x 17) satisfying (j)^ — on (0, T) x dfl, 



Q 



n 



U'<t't - f{u)[u^x - I{'^)]<l^xx - f'{u)ux{uxx - I{'^))<t>x dtdx 

UQ{x)4>(0,x)dx = Q. (10) 
Moreover, the function u satisfies for every t & [0, T], 

u{t,x)dx — / uo{x)dx, 



£Ht))+ I I f{u)[iuxx~Iiu))x]^dsdx < Siuo), (11) 
Jo Jn 

G{u{t))dx+ / / (uxx)'^ + Uxl{u)x dsdx < / G{uo)dx. (12) 
fi Jo Jn Jn 



We point out that the weak formulation (|T0)) involve two integrations by 
parts. Our second main result is concerned with non-negative initial data whose 
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entropies are possibly infinite (this is the case if uq vanishes on an open subset 
of O and n > 2). In that case, only one integration by parts is possible, and the 
solutions that we construct are weaker than those constructed in Theorem [TJ 
In particular, the equation is only satisfied on the positivity set of the solution 
and the boundary conditions are satisfied in a weaker sense. 

Theorem 2. Assume n > 1 and let uq G H^{fl) be such that uq > 0. For all 

T > there exists a function u{t,x) > such that 

ueC{o,T;L^{n))r\L°°{Q,T;H\n))r\C^'^n X (o,r)) 

such that 

and such that, for all (j) € l?([0,r) x Cl), 

II u(j)tdtdx+ // f{u)[uxx-I{u)]x(l)xdtdx+ / ua{x)cl>{0, x) dx = (13) 
JjQ JJp Jn 

where P = {{x,t) Cz Q : u{x,t) > {),t > 0}. Moreover, the function u satisfies 
the conservation of mass and the energy inequality 

Finally, vanishes at all points {x,t) of dCl x (0,T) such that u{x,t) ^ 0. 

Comments. These results arc comparable to those of U when A = 0. The 
reader might be surprised that they are presented in a different order than in 
[4]. The reason has to do with the proofs; indeed, in contrast with [4], weak 
solutions given by Theorem [2] are constructed as limits of the solutions given by 
Theorem[T] This is because the entropy is needed in order to construct the non- 
negative solutions. See the discussion at the beginning of Section [S] for further 
details. 

As pointed out earlier, the non-linearities in front of the stabilizing {uxxx) 
and destabilizing {{I(u))x) are the same. This is in contrast with the work of 
Bertozzi and Pugh [H [8] . By analogy with ([3]) , one could consider the equation 

Ut + [u^'Uxxx - U"\l{u))x)x = 

in which case a scaling analysis similar to that of [7] suggests that blow up 
can only occur if m > n -I- 1. However, to our knowledge, there is no physical 
motivation for such a generalization (in our case). 

Organization of the article. In Section [21 we give more details about the 
physical model leading to (|4]). We gather, in Section[3l material that will be used 
throughout the paper. In particular, we detail the functional analysis framework 
and the definition of the non-local operator / (which is similar to that used in 
[15]). Section m [5] and |6] are devoted to the proofs of the main results. Finally, 
we give in Appendix a technical result which is more or less classical. 
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2 Physical model 

In this section, we briefly recall the derivation of (jlj (see [19] for further de- 
tails). We consider a viscous liquid film which completely wets a solid horizontal 
substrate and is constrained between two solid walls (at a; = and a; = 1), see 
Figure 1. The fluid is Newtonian and is assumed to be a perfect conductor. The 




Figure 1: A viscous thin film submitted to an electric field E and gravity g 



substrate is a grounded electrode held at zero voltage. Thanks to the presence 
of another electrode (at infinity) , an electric field E is created which is constant 
at infinity (in the direction perpendicular to the substrate): 

E(a;,?;) — ^(0,£'o) as y ^ +oo. 

The height of the fluid is denoted by u{t,x). Under the assumptions of the 
lubrication approximation, it is classical that the evolution of u is described by 
Poiseuille's law: ^ 

ut - (jjl^^pj (14) 

where p is the pressure at the free surface of the fluid y ~ uit, x). This pressure 
is the sum of three terms: 

1. The capillary pressure due to surface tension, which can be approximated 

by 

(replacing the mean curvature operator by the Laplacian). 

2. The effect of gravity, given by 

P2 = gu. 

3. The additional pressure due to the action of the electric field E. 
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To compute the third term appearing in the pressure, we introduce the potential 
V such that E = —W, which satisfies 

AV = for y > u{x) 

and 

V{x, y) = on y = u{x). 
The condition at y — > c» means that wc can write 

V ^ EoiYo - y) 

with (using standard hncar approximation) 

AFo = for y > 0, .T e 

VFo ^0 asy-^oo, xefi (15) 

Yo{x, 0) = u[x) for X &VL. 

At the boundary of the cyhnder, we assume that the electric field has no hori- 
zontal component: 

dxV = 0, for X e 917, y > 0. 
The pressure exerted by the electric field is then proportional to 

P3 = 7Ey = ~-idyV{x,{)) = -lEoidyY^ - 1). 

The application u i— >■ 9^10(2;, 0) is a Dirichlet-to-Neumann map for the har- 
monic extension problem ^5]). We denote this operator by I{u). We will see 
in Section [3] that I{u) is in fact the square root of the Laplace operator on the 
interval Q, with homogeneous Neumann boundary conditions. 
We thus have 

P = Pi + J52 + = -lyUxx + gU~ jEoIiu) + Co 

for some constant cq, and we obtain ([1]) with c = f^, a = and A = 
Note that Poiseuille's law p4)) is obtained under the no-slip condition for the 
fluid along the solid support. Other conditions, such as the Navicr slip condition 
leads to 

f{u) = + Au-' 

with s = 1 or s = 2. This explains the interest of the community for general 
diffusion coefficient f{u). 

Boundary conditions. Along the boundary dft, the fluid is in contact with 
a solid wall. It is thus natural to consider a contact angle condition at x = 
and X = 1: Assuming that the contact angle is equal to 7r/2, we then get the 
boundary condition 

Ux — on 90. 
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In [in], the authors derive their analytic results in a periodic setting which 
is obtained by considering the even extension of u to the interval (—1,1) (recall 
that J7 = (0, 1)) and then taking the periodic extension (with period 2) to M. 

Finally, since the equation is of order 4, we need an additional boundary 
condition. We thus assume that u satisfies the following null-flux condition 

u^{uxx ~ I{u))x — on dil 
which will guarantee the conservation of mass. 



3 Preliminaries 

In this section, we recall how the operator / is defined (see [15]) and give the 
functional analysis results that wc will need to prove the main theorem. A 
very similar operator, with Dirichlct boundary conditions rather than Neumann 
boundary conditions, was studied by Cabrc and Tan [5]. 

3.1 Functional spaces 

The space Hf^{ft). Wc denote by {A^-, (pk}k=Q.i.2... the eigenvalues and corre- 
sponding cigcnfunctions of the Laplace operator in 51 with Neumann boundary 
conditions on dfl: 

-Aifk = Xky^k in . . 

d,^ipk =0 on dfl, ^ ' 

normalized so that (fj. dx = 1. When Q, = (0, 1), we have 

Ao = 0, (po{x) = 1 

and 

Afc = (fcTr)^ , <y3fc(x) = \/2cos(fc7ra;) fc = 1, 2, 3, . . . 

The (/Pfc's clearly form an orthonormal basis of L^iVl). Furthermore, the ipkS 
also form an orthogonal basis of the space Hf^{n) defined by 



H-^{n) = lu = J2 m ; ^ c2(l + A^.) < +(» \ 

{ k=0 k=0 ) 

equipped with the norm 

oo 
k=0 

or equivalently (noting that cq = u{x) dx and A^ > 1 for k > 1): 

= l|w!lL(fi) + !l"ll|^(o) 
where the homogeneous norm is given by: 



fc=i 
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A characterisation of iJ|^(r2). The precise description of the space Hf^{n) 
is a classical problem. 

Intuitively, for s < 3/2, the boundary condition = docs not make sense, 
and one can show that (sec Agranovich and Amosov ]T] and references therein): 

H'^{n)^H'{n) forallO<s<^. 

In particular, we have H^{V,) = H^{n) and we will see later that 

/ / {u{y)- u{x)fv{x,y)dxdy 
ti^y^^ JnJn 

where v{x^ y) is a given positive function; see ()20p below. 

For s > 3/2, the Neumann condition has to be taken into account, and we 
have in particular 

H%{n) = {ue H^{VL) ; = on dVt} 

which will play a particular role in the sequel. More generally, a similar charac- 
terization holds for 3/2 < s < 7/2. For s > 7/2, additional boundary conditions 
would have to be taken into account, but we will not use such spaces in this 
paper. In Sectional we will also work with the space Hff{n) which is exactly 
the set of functions in H^{il) satisfying = on dfl. 

The case s = 3/2 is critical (note that u,^\gn is not well defined in that space) 
and one can show that 

where d{x) denotes the distance to dfl. A similar result appears in [3]; more 

3 

precisely, such a characterization of H^(rt) can be obtained by considering 
functions u such that G Vo(ri) where Vo(f^) is defined in [9] as the equivalent 

1/2 

of our space (il) with Dirichlet rather than Neumann boundary conditions. 
We do not dwell on this issue since we will not need this result in this paper. 

3.2 The operator / 

As it is explained in the Introduction, the operator / is related to the compu- 
tation of the pressure as a function of the height of the fluid. 

Spectral definition. With Xk and (pk defined by (|16p , we define the operator 

oo oo ^ 

/:^Cfe(^fei — > -"^CkX^ipk (17) 

k=0 k=0 

which clearly maps H^{^) onto L^(ri) and Hff{n) onto H^{^). 
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Dirichlet-to-Neuman map. We now check that this definition of the op- 
erator / is the same as the one given in Section [2l namely / is the Dirichlet- 
to-Neumann operator associated with the Laplace operator supplemented with 
Neumann boundary conditions: 

We consider the following extension problem: 

-Aw = infix (0, +00), 

v{x,0) u{x) onQ, (18) 
Vi, = on dfl X (0, 00). 

Then, we can show (see [U]): 

1 

Proposition 1 (|15j). For all u G H^{fl), there exists a unique extension 
V e H^{nx (0, +00)) solution of J7?)j. 

Furthermore, if u{x) = 'YlV=i^kVh{x), then 

00 

. ^ 1 

v{x,y) = 2_,CkVk{x)c-yi-p{~\ly). (19) 

k=l 

and we have: 

Proposition 2 ([15]). For all u E Hjf{n), we have 

dv 

I{u){x) = — — (.T,0) = dyv{x,0) for all x Eft, 

where v is the unique harmonic extension solution of I118\) . 
Furthermore I o I[u) = —Au. 

Integral representation. Finally, the operator / can also be represented as 
a singular integral operator: 

Proposition 3 (|15j). Consider a smooth function u : f2 — > M. Then for all 

X efi, 

I{u){x)= / {u{y) -u{x))v{x,y)dy 



where u{x,y) is defined as follows: for all x,y d fl, 

^^^2 (l-cos(^(x-y)) + l-cos{n{x + y))) ' 
3.3 Functional equalities and inequalities 

Equalities. The semi-norms || • j^,, || ■ II ' ILi,„, ^^d || • ||^f2 (^j) 

are related to the operator / by equalities which will be used repeatedly. 

Proposition 4 (The operator / and several semi- norms - |15j). 
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For all u G H 2 (n), we have 



-J ul{u)dx = \ j^jj^^^^) ~ u{y)fv(x,y)dxdy = ll",,^^^^) 



|2 



where v is defined in (j20p . 
For all u S Hj^iVl), we have 



Uxl{u)x dx = llwlP 3 



For all fc e N and u G H^+'{n), we have 

^{d^I{u)fdx=\\u\\l,^^,^^y (21) 

Inequalities. First, we recall the following Nash inequality: 

It implies in particular that, 

MlHn)<l\Mhin)+CMlHny (22) 

This inequality will allow us to control the norm by the energy £{u) and 
the norm. Indeed, we recall that the energy is defined by 

= llj^.f+ul{^)dx^^\\u\\l,^^^^ - l\\u\\li^^^y (23) 

We then have: 

Lemma 1. There exist positive constants a, /3 such that for all u G H^{i}), 

Remark 1. See also Lemma 4.1 in |19| . 
Proof. We have 

\\u\\l,,^^^2£{u) + \\u\\^. , , 
II iiifi(o) \ ' II "H^{n) 

and using (PT|) with fc = and (P^. we get: 

lliill^. 1 = — u I(u) dx 



< 



\u\\LHn}\\Iiu)\\mn) 



< ll"l|L2(n)||w||_ffi(o) 



< 



^lklli^(o) + ^lkll|i(n) 

hence the result. □ 
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4 A regularized problem 

Wc can now turn to the proof of Theorem [1] and [H As usual, we introduce the 
following regularized equation: 



which satisfies 

£<fe{u)<M foraUuGR. 
Ultimately, we will show that the solution u satisfies 

< u{t, x) < Mo 

for some constant Mq independent of A/, so that wc do not have to worry about 
M (provided we take it large enough). The e is of course the most important 
parameter in the regularization since it makes (|24p non-degenerate. 

The proof of Theorem [1] consists of two parts: First, we have to show that 
the regularized equation ([M)) has a solution (which may take negative values). 
Then we must pass to the limit e — > and show that we obtain a non-negative 
solution of dl]). 

In this section, we prove the first part. Namely, we prove (to be compared 
with Theorem 1.1 in [4]): 

Theorem 3. Let uq E H^{n). For all T > there exists a function u'^{t,x) 
such that 




i{u))x)x = for X en, t>0 

- I{u))x =0 for xedn, t>0 



(24) 



where the mobility coefficient f{u) is approximated by 



fe{u) 



min(max(e,/(|u|)),M) 



such that, for all (j) e I?([0,T) x n), 




(25) 



Moreover, the function satisfies for every t € [0,r], 




(26) 
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and 



{ulfdx+ f [ Uu-){ul^,fdsdx 
JO Jn 

■t 

,2 



< / {{uo),rdx+ I / f,iu')ul,,{Iiu')),dsdx (27) 
Jn 



and 



t 

£ \2 



Ge{u^t))dx+ / / {ul^y +ull{u^)^dsdx < / Geiuo)dx. (28) 
^0 Jn Jn 

where is a non-negative function such that /gG" = 1. 

Finally, is ^-Holder continuous with respect to x and ^-Holder continuous 
with respect to t; more precisely, there exists a constant Go only depending on 
ft and \\u^\\L^(Q,T;m(n)) and ||/e(u^)[u^^ ~ I{u^)]x\\l^(q) such that 

M\\ 11 < Co. (29) 

Remark 2. This theorem is very similar to Theorem 1.1 in [4], and some steps in 
our proof follow along the lines of [4]. For instance, getting the Holder estimates 
from the L°°H^ estimate is done in the same way. However, the main difficulty 
in proving this theorem is precisely to get the L°°H^ estimate; this step is not 
straightforward at all and this is a significant difference with [4]. 

Proof of Theorem\^ Theorem [3] follows from a fixed point argument: For some 
T*, we denote 

V = L'{Q,T,-H%{fl)) 

and we define the application F : V ^ V such that for v ^ V , J-{v) is the 
solution u of 

Ut + {fe{v){Uxxx - I{v)x))x =0 iOTXeil, t>0 

ux = 0, fe{v){uxx-l{v))x=0 ioTxedn, t>Q (30) 
u{x, 0) ~ uo{x) for X g 57. 

The fact that T is well defined follows from the observation that for v £ V, 
we have 

ait, x) = fe{v{t, x)) € [e, M] and g{t, x) = I{v)x € L^{Q) 
and the following proposition: 

Proposition 5. Consider uo E H^{n) and a{t,x) G L°°{Q) such that e < 
a(t,x) < M a.e. in Q. If g G L'^{Q), then there exists a function 

u G C(0, T; L^{VL)) n i°°(0, T; H^{n)) n L^{Q, T; H%{n)) 
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and 

\u(f>t + a{uxxx — g)4'x\ dt dx + / uo{x)<j){0, x) dx = 



for all (j) € I?([0, T) x Q). Moreover, for every t G [0,T], u satisfies 

u{t,x)dx= / uo{x) dx 
and 

{ux f{t)dx + l- ( I a{uxxxY ds dx < ^ ( ( g^dsdx+ / (uq)^ da;. (31) 
^ Jo Jn ^ Jo Jn Jn 

Furthermore, u is ^-Holder continuous with respect to x and ^-Holder continu- 
ous with respect to t; more precisely, there exists a constant Co only depending 
on n and |iM|lL~(o,T;_Hi(s^)) \Wxxx - g\\L^(Q) such that 

\\u\\ 11 <Co. (32) 

This proposition is a very natural existence result for the fourth order linear 
parabolic equation 

ut + {auxxx)x = iag)x- 
Its proof is fairly classical, wc give some details in Appendix |X] for the interested 
reader. 



Next, we show the following result: 

Lemma 2. There exists a (small) time > 0, depending only on e, M and 
n, such that T has a fixed point u in V = L^{0, T^,; Hfj(n)) for any initial data 
Uo G H^{Q,). Furthermore, u satisfies 

\\u\\v < R\\uo\\Hi{n) 

and 

ll"llL-(o,T.;_fri(n)) < ^/2||uo||^i(f2). (33) 
Before proving this lemma, let us complete the proof of Theorem |3l 

Construction of a solution for large times. Lemma [2] gives the existence 
of a solution u\ of ([24]) defined for t G [0, T,]. Since T, does not depend on the 
initial condition, we can apply Lemma [2] to construct a solution in [T*,2r*] 
with initial condition u\(T^,x) which is H^ijl) by (|33p . This way, we obtain a 
solution of ([M)) on the time interval [0, 2T,]. Note that we also have 

!l"''llL°°(o,2T.;ffi(n)) - ^ ll"o||/i-i(n)- 

Iterating this argument, we construct a solution on any interval [0,r] satis- 
fying in particular, for all /c G N such that kT^ < T, 

ll"''llL°°(0,fcT.,ffi(n)) - ^ -R|l'"o||_H-i(n)- 
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Energy and entropy estimates. The conservation of mass follows from 
Proposition [SJ but we need to explain how to derive ((26)) . ((27| and p8|) from 
(|^ . Formally, one has to choose successively (j) = — + = —u%^ and 

= Gg(u'^). Making such a formal computation rigourous is quite standard; 
details are given in Appendix [B] for the reader's convenience. Finally, (|29p 
follows from □ 

Proof of Lemma\^ We need to check that the conditions of Leray-Schauder's 
fixed point theorem are satisfied: 

F is compact. Let be a bounded sequence in V and let w„ denote 

F{vn). The sequence {I{vn))x is bounded in L'^{Q), and so 

5n = fe{vn)dxl{vn) is boundcd in L^iQ). 

Estimate (|3T|) implies that u„ is bounded in i^(0, T*; In particular 

dxxxUn is bounded in L^(0, T*; L^(f2)) and Equation ([50)1 implies that dt{un) is 
bounded in L^(0, T*, H~^(r2)). Using Aubin's lemma, we deduce that (u„)„ is 
pre-compact in F = L'^{Q, T^;H%{n)). 

T is continuous. Consider now a sequence in V such that w„ — > w in 

y and let u„ = T{vn)- We have in particular in L^{Q) and, up to 

a subsequence, we can assume that u„ — > u almost everywhere in Q. Hence, 
fe{vn) fe{v) almost everywhere in Q. We also have that {I{vn))x converges 
to {I{v))x in L'^{Q), and since \fe{vn)\ < M a.c., we can show that 

9n = fe{Vn)dxI{Vn) fe{v)dxl{v) = Q in {Q) . 

Next, the compacity of T implies that is pre-compact in the space 

L'^{0,T; Hjf{V,)), and so u„ converges (up to a subsequence) to U in V. In 
particular, u„ — J7 in L^{Q) and (up to a another subsequence), w„ — > [/ 
almost everywhere in Q. We thus have fe{un) fe{U) in L^{Q), and passing 
to the limit in the equation, we conclude that U = u = J-{v) (by the uniqueness 
result in Proposition [5]). Since this holds for any subsequence of u„, we deduce 
that the whole sequence u„ converges to u hence 

J^(vn) J^{v) in y as n — oo 

and J- is continuous. 

A priori estimates. It only remains to show that there exists a constant 
R > such that for all functions u G V and a G [0, 1] such that u = aF{u), we 
have 

\\u\\v < R- 

This is where the smallness of will be needed. 
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Using ([5T|). we see that 



o 2 jQ 

<^/ niliu),j'dxdt+ Huojidx, (34) 
and using (j2ip and the interpolation inequality 

||":r|lL2(o) < C\\u\\l2(^^-j\\Ua;x\\l2(^i^y (35) 

we get: 

' / i{I{u)^fdxdt < 2C, [ ' \\Iiumf2dt+^ [ ' \\iliu)),,{t)\\ldt 
Jn Jo Jo 

< 2C, / ||M,(t)||2di+— / \\u^^4t)\\ldt 

i-T, 



2Af JO 



< 2C,r* sup \\u,{t)\\i + — ||Mx.x(t)||2dt(36) 
te[o,T,] 2A:; Jo 



where only depends on the constant C in ((35|) and the parameters e and M. 
Combining (IMl) and we conclude that 



(l-MCeT*) sup ||M,(t)||^ + ^ / * ||u,,,(i)|l2dt < / iuo)ldx. 
te[o,T] 4 Jo Jn 

Therefore, choosing JMTT' S*^^ following estimates 

2 

ll"llL°'(0,T.;ffi(a)) - ^ll'"o|liji(n) a-nd ll'^!lL2(0,T.;ff3,(O)) - ^ 1 1^0 1 li^i (n) 



Since we also have 



u{t) dx = j uo dx, 
n Jn 



wc deduce that \\u\\v < R for some constant R depending on e, which completes 
the proof. □ 

5 Proof of Theorem [1] 

As pointed out in the introduction, one of the main difficulties in the proof of 
Theorem[T]is that the natural energy estimate ([5T]) docs not give any information 
by itself, since S{u) may be negative. Even if Lemma [T] implies that the quantity 
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is bounded below by the norm of u, the mass conservation only allows us to 
control the norm of if we know that u'^ is non-negative. Unfortunately, it 
is well known that equation (|24p does not satisfy the maximum principle, and 
that the existence of non-negative solutions of Q is precisely a consequence 
of the degeneracy of the diffusion coefficient, so that while we can hope (and 
we will prove) to have lim^^^o '>^^ > wc do not have, in general, that u"^ > 0. 
Lemma [3] below will show that it is nevertheless possible to derive some a priori 
estimates that are enough to pass to the limit, provided the initial entropy is 
finite (O. 

Proof of Theorem [H Consider the solution of (j24p given by Theorem [31 In 
order to prove Theorem [U wc need to show that lim5_j.o u"^ exists and solves 

m- 

Since wc cannot use the energy inequality to get the necessary estimates on 
u^, we will need the following lemma: 

Lemma 3. Let denote the following functional: 

He{v)= I [vl + 2MGe{v)]dx. 
Jn 

Then the solution given by Theorem satisfies 

H,{u%t)) + ^ f I {u%.,fdxds + \ f I f,{un{ul,J'dxds<HMe'/' 
Jo Jn ^ Jo Jn 

for every t G [0, T] . 

Proof of Lemma[M Using ([27]) and (f28| we see that 

/ G,{u^{t))dx + l f I [ulJdxds 
Jn ^ Jo Jn 

< f G,{uo)dx + l- f j {ulfdxds 
Jn ^ Jo Jn 

and 

/ /* / f,{u-){ul^^fdxds 

Jn ^ Jo Jn 

< I {{uo)^fdx + ^ f f {ul^fdxds. 
Jn ^ Jo Jn 

This implies 

He{u'{t))<He{m)+ I He{u'{s))ds, 

Jo 

and Gronwall's lemma yields the desired result. □ 
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Sobolev and Holder bounds. We now gather all the a priori estimates: 
Using the conservation of mass, Lemma [3] and inequahty (|9]), we see that there 
exists a constant C independent of e such that: 



sup / [ul{t)fdx < C, (37) 

t<£[0,T] Jn 

sup / G,{u^{t))dx < C, (38) 
r [ Mu')[u%,jUxdt < C, (39) 

f f {uU'dxdt < C. (40) 
Jo Jn 

Next, we note that (|37p yields 

£{u^) > -||"''llL=o(o,T;iji/2(o)) > -C||u^||Loo(o,T;_f/i(n)) > -C 
and so gives 

/ / Mu%ul,,-Iiun,,]'dsdx<C. (41) 
Jo Jn 

Finally, estimates ^ and ^ yield that u" is bounded in C^^f '^/^(Q). 

Limit £ — 0. The previous Holder estimate implies that there exists a func- 
tion u{x,t) such that converges uniformly to u as e goes to zero (up to a 
subsequence). Inequality (|40|) also implies that 

^uin L^{0,T;H^{n))-wcak 

and Aubin's lemma gives 

u'' — >u in L^{0,T;H\n))-stTong. 
After integration by parts, (|25p can be written as 



u-</>i - A(u^)[uL - I{u-M., - au')ul{ul, - I{u-))(p, dtdx 

I 

= / uq{x)4>{Q,x) dx 
Jn 

and passing to the limit e — >■ gives pI7|) . 

Non-negative solution. It only remains to show that u is non-negative. This 
can be done as in [4], using ((38)) (and the fact that / satisfies ([5]) with n > 1). 
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L°° a priori estimate. Finally, ([57)) and Sobolev's embedding implies that 
there exits a constant A/q depending only on ||ito||ffi(f2) such that 

Q<u{t,x)<Mo. (42) 

Choosing M > AIq, we deduce that u solves (|3]). □ 

6 Proof of Theorem [2] 

In order to get Theorem [21 we need to derive the following corollary from The- 
orem [TJ 

Corollary 1. The solution u constructed in Theorem [1\ satisfies for all 4> G 

p((o,r)xf]), 

u(j)t dtdx+ / / /(u) [w^rx - I{u)]x4>x dtdx = (43) 



p 

where P = {{x,t) G Q : u{x,t) >0,t> 0}. 

Proof. In view of the proof of Theorem[Tl u is the uniform limit of a subsequence 
of (u^)e>o where is given by Theorem [31 Since u^ satisfies ([^5]) . it is thus 
enough to pass to the limit in this weak formulation as e — > in order to get 
the desired result. Let ft,^ denote fe{u^)[u'^.^ — I{u)]x. Estimates ([41]) and (|42p 
imply 

hi dx dt < C. (44) 

'Q 

In other words, {h'^)e is bounded in L^{Q). Hence, up to a subsequence, 

h'^ h in L^{Q)—wcak. 

Furthermore, we recall that there exists a continuous function u{x, t) such that 
converges uniformly to u as e goes to zero (up to a subsequence). 
Passing to the limit in ([25[) . we deduce that the function u satisfies 

u(j)t dt dx + j j h(j)x dt dx = 0. 



Wc now have to show that 



in {u = 0}, 

fiu)[u.,x-Iiu)]^ mP^{u>0}. 



First we note that for any test function cf) and > 0, we have 



fe{u')[ul.xx - Iiu'')x](j}dx dt 

J {«<!)} 



<C{cP)(f,{3v/2)){ [ [ f,{u')[ux,x~I{u),]''dxdt 



1/2 
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for £ small enough (so that |it^ — u| < 77/2). Inequality (PT|) thus implies 

< C(0)/(r;/2)i/2. 

We deduce (since /(O) = 0) 



lim sup 



T 

fe{u'')[ul^^ ~ I{u^)^]<l>dx dt 

{u<2ri} 



h^O on {w = 0}. (45) 
Next, pi)) yields (for e > and 77 small enough) 

\ul,,~I{u'),\^dxds<C{Ti). 

{li>2?7} 

This implies that, if Q,, denotes {u > 2r]}, (u^^.^, — I{u'^)x) is bounded in the 
space L^{Qn)- Hence, we can extract from (wj^.^, — I{u^)x)e>o a subsequence 
converging weakly in L^((5^). Moreover, remark that is an open subset of 
Q (recall that u is Holder continuous) and u'^^xx ~ Iiu'^)x converges in the sense 
of distributions to Uxxx — I{u)x (use the integral representation for /(•)). We 
thus conclude that, 

^Ixx - I{u^)x Uxxx - I{u)x in L'^iQri)- 

This yields 

h = f{u) [uxxx - I{u)xx] in {u > 0} 
which concludes the proof of Corollary [TJ □ 

Wc now turn to the proof of Theorem [2l 

Proof of Theorem\E When uq does not satisfy (0), we lose the L'^{0, T; H'^{n)) 
bound on u", and the previous analysis fails. However, we can introduce 

Uq = Uq + S 

which satisfies ([9]). Theorem [1] then provides the existence of a non-negative 
solution of (|10p . In view of Corollary [TJ satisfies: 



5, 

u 

Q 



u^<j>tdtdx+ // f{u^)[ui^~ I{u^)]x(t)xdtdx = 0. (46) 



Since is non-negative, the conservation of mass gives a bound in the space 
L°°{0,T; L^{n)) and allows us to make use of the energy inequality: Indeed, 
using and Lemma [T] we see that there exists a constant C independent of 
6 such that 



l"'^l|L°°(0,T;_f/i(n)) < C 



and 



I f f{u')[ui^x-Iiu%]'dsdx<C. (47) 
Jo Jn 
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We now define the flux 

^f{u%ui,,-I{u%]. 

Inequality (|47|) implies that is bounded in L^{Q), and that there exists a 
function h G L^{Q) such that 

in L2(Q)-weak. 

Proceeding as in the proof of Theorem [1] we deduce that is bounded in 
(ji/24/8^j^x (0, T)) and that there exists a function i) such that converges 
uniformly to u as ^ goes to zero (up to a subsequence). We can now argue (with 
minor changes) as in the proof of Corollary [1] and conclude. □ 

A Proof of Proposition [5] 

Our goal here is to prove the existence of a weak solution of 

Ut + {auxxx)x = (ag)x- 
We flrst prove the following proposition. 

Proposition 6. For all h E H^{n), there exists v G Vq := fl such that 
for all 4> e V{Vl), 

V — h 



dx + / avxxx4>x dx = ag<j>x dx. (48) 



n 'T J Jn 



V dx ~ h dx, 
n Jn 

^ ^ 2 , _ /" „„,2 ^ 1 / ;,2 



In particular, 



V-+T I av-^.-^ < - h^ + T agvxxx- (49) 
n Jn ^ Jn Jn 

Proof. In order to prove this proposition, we have to reformulate the equation. 
More precisely, instead of choosing test functions (p e 2?(f2), we choose (p = 
—i^xx + dx where ip is given by the following lemma 

Lemma 4. For all cj) S ^{^1), there exists ip £ ^{Vl) such that 

-i^xx + / ipdx ^ (p. 
Jn 

Hence, wc consider Vq := f] Hf^ equipped with the norm HwHyj, = 
||fa;2-a;||/,2 + (/ dx)"^ and wc look for ?j G Vo such that for all -ip G 'D{Cl), 

Vxtpx dx + T avxxxipxxx dx + i vdxj i i'dx 
n J \Jn / \Jn 

= / hxipx dx + i / hdxj i / ipdxj +T agtpxxx dx. (50) 
Jn \Jn J \Jn J Jn 
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We thus consider the bilinear form A in Vq defined as follows: for all w, u; S Vq, 



A{v,'w) = / VooWxdx + T / aVxxxWxxx dx + / vdx\ I / wdx 
Jn J \Jn / \Jn 

We check that it is continuous and coercive: 

< ||Wa;||2||W:r||2 + A/r||Uj;j;j;|l2||Wxxa;||2 + 

< C||i;||yJ|w||yo 



A{v,v) > [ivxr + eT{vxxxy]dx+[ vdx] > \\v\\t,^. 
Jn \Jn J 

We now consider the following linear form L in Vq: for all w e Vbi 



L(w) = I hxWx dx + / /ida:) / w dx ] + t agWxxxdx 
Jn \Jn J \Jn J Jn 

Since < a < M and g E L^, L is continuous as soon as h E H^{fl). Lax- 
Milgram theorem thus implies that there exists v E Vq such that (|50l) holds true 
for all w EVq. 

Eventually, remark that conservation of mass and are direct conse- 
quences of pS)) . The proof of Proposition IH] is now complete. □ 

We can now prove Proposition [5l 

Proof of Proposition^^ For any r > 0, we consider Nj- = [t"]- We then define 

inductively a sequence (u")„=o of Vq as follows: m° = uq and u"'^^ is 

obtained by applying Proposition[6]to h — u". We then define u'^ : [0,NrT) x ^l 
as follows: 

u''{t,x) = u"-{x) for t e [nr, (n + l)r). 

We have u'^{t, x) dx = uo{x) dx for all t. We also derive from (|49p that we 
have 

{ulf{T,x)dx+ / / a{ulxxf{t,x)dtdx 
Jo Jn 

< / {{uo)x)'^dx+ / / {aguxxx){t,x)dtdx 
Jn Jo Jn 

In particular, {u'^)r is bounded in L°°(0, T; iJ^^)) and (5^ — u^)t is bounded 
in L2(o,T - r;7J-i(f^)) where Srv{t,x) = v(t + t,x). We derive from [HI 
Theorem 5] that (m^)t is relatively compact in C(0, T; 

We now have to pass to the limit in (gS]). Since {ul.^^)r is bounded in L'^{Q) 
and we can find a sequence t„ — such that u^^ u in C(0, T, i^(rj)) and 
''^''xxx ~^ Uxxx in L^{Q). This is enough to conclude. 

We next explain how to get ([32]) . Sobolev's embedding imply that there 
exists a constant K (depending on \\u\\i:^<xi(^Q X;H^{n))) such that 

\u{xi,t)~u{x2,t)\ < K\XI-X2\^'^ 
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for all xi, X2 £ ^ and for all t E (0,T). Since u satisfies 

Ut = K 

with h G L^{Q), it is a fairly classical result that Holder regularity in space 
implies Holder regularity in time. More precisely, we have (see [J, Lemma 2.1 
for details): 

Lemma 5. There exists a constant C such that for all x\, X2 in f2 and all ti, 
t2 > 0, 

\uixut,) - Uix2,t2)\ < C\XI - X2\^l^ + C|ti - fsT/'. 

The proof of Proposition [S] is thus complete. □ 



B Proof of ([2SD, (I27D and (EHD 

We have to derive (UHl), (EB and (US]) from (gS]). Using the fact that u"" lies in 
C([0, T], iJ^(i7)), we first state the following lemma 

Lemma 6. For all (p € 'D{Q), 

jj u^(t)t + feiu'')[ul^~ I{u^)]x(j)xdtdx 

= I u^{x,T)(j){x,T)~ [ uo{x)(f>{x,0)dx. (51) 
Jn Jn 

The proof of such a lemma is fairly classical. It relies on mollifiers that are 
decentered in the time variable. More precisely, one considers a smooth even 
function p : R — > [0, 1] compactly supported in [—1,1] and such that / p = 1. 
Then for a > and 5 G M, one can define 

Pa,s{t) = Tspa{t) = ap 

If now a function / is defined in [0,r], it can be extended by to R; in other 
words, it can be replaced with /Iq where la{x) = 1 if x S and Inix) = if 
not; then the convolution product in R: {fin) * Pa.s is a smooth function in R 
which vanishes near t ^ (resp. t ^ T) ii 6 > (rcsp. (5 < 0). 

Consider a smooth function p such as in the proof of Lemma [SI Consider 
a > and define pa{x) = ap{^). Then consider 0{x,t) = pa{x)pa(t). 

Lemma 7. Recall that u" e L°°iQ) and h" = /e(u=)[u^^ - I{u^)]x e L'^iQ)- 
Then for all v G L^{Q), 

JJ u'iv*e)t^ JJ iu'i.9)tv (52) 

where f * g means (fig) * (gig). 

We next apply Lemma [6] with (f) = v -kO where v is chosen to be successively 
u^^ -kd, -kO) and G'^{u^ -kO). After direct computations, we can let a — > 
and get the desired estimates. 
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